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t-H Abstract. We analyze the features of the Minkowskian limit of a particular non- 

analytical f(R) model, whose Taylor expansion in the weak field limit does not hold, 
as far as gravitational waves (GWs) are concerned. We solve the corresponding Einstein 
equations and we find an explicit expression of the modified GWs as the sum of two 
f"**) terms, i.e. the standard one and a modified part. As a result, GWs in this model are 

t-H not transverse, and their polarization is different from that of General Relativity. The 

velocity of the GW modified part depends crucially on the parameters characterizing 
the model, and it mostly results much smaller than the speed of light. Moreover, this 
investigation allows one to further test the viability of this particular f(R) gravity 
?_i theory as far as interferometric observations of GWs are concerned. 



PACS numbers: 04.50.Kd Modified theories of gravity; 04.30.-w Gravitational waves; 
04.25.Nx Post-Newtonian approximation, perturbation theory, related approximations; 
04.20.Jb Exact solutions; 04.20.Cv Fundamental problems and general formalism 
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1. Introduction 

Within the possible modifications of General Relativity (GR), /(R) modified theories of 
gravity are based on replacing the Ricci scalar R in the Einstein-Hilbert (EH) action by 
a suitable function /(R) of it. Even though the predictions of GR reproduce extremely 
successfully the great majority of the phenomena within a very stringent accuracy, there 
are nevertheless some experimental data which are not matched by GR. Some of these 
phenomena can be explained as due to the presence of extra (otherwise unobserved) 
matter or energy contributions. /(R) theories of gravity are a tool to describe these 
phenomena as due to effects of the geometry (i.e. 'on the left-hand side of Einstein 
equations') rather than of some kind of matter (i.e. 'on the right-hand side of Einstein 
equations'). There are at the moment many analyses aimed at discriminating from an 
experimental point of view between geometrical effects and matter effects. 
It is then important to constrain the parameter space of these models in such a way 
to reproduce these unexpected phenomena without discarding all the other well-tested 
predictions of GR. Moreover, it is important to stress that the suitably-constrained 
parameter space should be valid at all physical scales. 

It is interesting to remark that some kind of /(R) theories of gravity descend from a 
suitable 'low energy limit' of some other unification theories. We remark that, in the 
latter case, the presence of other invariants (built as suitable products of the Riemann 
and the Ricci tensor as well) is also possible pQ. 

In this paper, we study the behaviour of the gravitational waves in the non-analytical 
/(R) model proposed in [2], i.e. /(R) = R + 7R^, since the corresponding weak 
field equations posses a peculiar feature, for which retaining non-linear terms in the 
dynamics makes sense. In fact, as far as the parameter 2 < < 3, there are non-linear 
corrections to the General Relativity analysis which, being greater than quadratic term 
and tending to the same order of the linear ones as /3 — > 2, can not be disregarded 
and have a significant impact on the theory predictions in the limit of weak spacetime 
ripples propagation. 

Our study outlines how, in parallel with the linear vacuum massless gravitational 
waves of General Relativity, a non-linear wave appears, having a non-trivial (non- 
transverse and non-traceless) morphology and whose amplitude increases with time. 
The effect induced by this "anomalous" wave on the test particle geodesic deviation is 
described in some detail, and the request that the modified tidal term be below the 
typical experimental constraints lead us to restrict the velocity range of these non-linear 
ripples. Indeed, a peculiarity of the considered model consists of the non-null character 
of the propagation vector, whose modulus results to be fixed by the initial and boundary 
conditions on the "radiation" field, and must take values corresponding to a propagation 
velocity well-below the speed of the light, apart from the limiting case — > 2. 

Finally, we provide a study of the polarization of the modified wave, based on the 
action of the tidal force on a particle system initially arranged on a circumference and 
having as reference the comparison to the standard case. This characterization of the 
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anomalous deformation of the polarization ellipses offers a valuable tool for recognizing 
the presence of this non-linear feature among the background signals detected by ground 
interferometers. In particular, the increasing power-law behaviour that the modified 
contribution outlines at a fixed point of space makes this new effect, whether existing, 
very promising for a future detection. In this sense, the present study, on the one hand, 
offers a very reliable trace to identify a specific gravitational dynamics modification, 
and, on the other hand, offers an intriguing scenario for enhancing the detection of 
gravitational signals, which, in the considered case, are particularly able to bring 
information on their generating source, here schematized by the initial and boundary 
condition on their propagation. 

The paper is organized as follows. In Section[2j we briefly recall the basic statements 
about /(R) theories of gravity, with particular attention to their weak- field limit. In 
Section |3j we review the main features of a particular class of non- analytical /(R) 
theories of gravity, i.e. /(R) = R + tR' 3 , where the standard EH action of GR is 
perturbed by the addition of a term, for which the Taylor series in the vicinity of R — 
does not hold. In section |4j we first briefly review the main features of the weak-field 
limit of /(R) models, analytically solve and discuss the weak-field limit of the Einstein 
equations for this model, and explicitly write down the form of gravitational waves; to 
do so, we analyze the physical meaning of the degrees of freedom of the model. We 
also compare our results with those obtained in the case of analytical /(R) models. In 
Section |5j we discuss the interaction of gravitational waves with test particles. Brief 
concluding remarks end the paper. 

2. /(R) modified theories of gravity 

/(R) theories of gravity are obtained when the Ricci scalar R in the Einstein-Hilbert 
action 

is replaced by an arbitrary function /(R) of it, such that the modified gravitational 
action reads 

S = J *xV=sf<&, (2) 

where y/— g is the determinant of the metric tensor g^, a/— g = detg^ u . In vacuum, 
the modified Einstein equations take the form 

-^/(R) + f (R)R M , - V„V„/'(R) + g^VpV.f (R) = 0, 

3g^V p V CT /'(R) + /'(R)R - 2/(R) = 0, 

where R^j, is the Ricci tensor, is the covariant derivative, and a prime denotes 
differentiation with respect to R, /'(R) = df/dR. Here and in the following, we will 
use bold letters for the full covariant objects (such as g M „, R^, R), and usual letters for 
their weak-field expression (such as g^, R, R pv ). We will adopt the following notation: 



Non analytical f(R) theory of gravity 



4 



Greek indices run form to 3, i.e. /j,, v — 0,1,2,3, Latin indices run form 1 to 3, i.e. 
i,j = 1, 2, 3, and we follow the standard notation [3], with mostly-minus signature. 

The physical effect of f(R) models is to add a scalar degree of freedom to EH gravity 
(see, for example, the review [I] and the references therein). In the Jordan frame, the 
dynamical features of the non-constant first derivative f'{R) can be treated as a scalar 
field. In the scalar-tensor version of the models, via a conformal transformation, the 
modified action ^ is rewritten in terms of a scalar field minimally coupled to gravity 
in the Einstein frame. 

General weak-field limit In the weak-field limit, in vacuum, we consider the metric 
tensor g^ u as consisting of the flat metric rj^ = (1, —1, —1, —1) and a small perturbation 
of it, \hfj, v \ « 1, i.e. 

9ixv = + h^u- (4) 

The weak-field limit of Einstein equations consists in considering the 0(/i)-terms only, 
and in neglecting higher-order (9(/i 2 )-terms. 

If the function /(R) admits a Taylor expansion in the neighborhood of R ~ 0, its 
Taylor series f{R) = X]j=o° a>jR? (where the a^s are the j-th order Taylor coefficients of 
the series and have the dimension of length 2 i~ 2 ) can be inserted in the modified Einstein 
equations. 

Even though an infinite number of parameters has, in principle, to be fixed, the weak- 
field limit of this class of /(R) models can be shown to depend only on the parameters a\ 
and Gt2, for vanishing a®. In fact, the limit of flat Minkowski space (in GR) is recovered 
for R ~ 0, such that ao ~ 0. 

Furthermore, we remark that, even though any value of the parameter a\ is in principle 
admitted from the Taylor expansion, if the Taylor series of f{R) has to be interpreted 
as a perturbation of the Ricci scalar in the standard EH action, then the value of a\ 
should be close to 1. In fact, the effect of a\ ^ 1 can be interpreted, in the comparison 
with GR, as a modification of the value of the gravitational constant G. (See [5], [6]). 

If, on the contrary, the function /(R) is not analytical around R ~ (i.e. its 
Taylor expansion does not hold in the neighborhood of R ~ 0), other paradigms are to 
be looked for. It is within this perspective that non-analytical /(R) models have been 
addressed. 

For example, the model /(R) = R — fi 2 sin (where fi is a constant with the 
dimensions of length^ 1 and A is a constant with the dimensions of length -2 ) has been 
proposed in [7]. 

Furthermore, the /(R) model consisting of a sum of different powers of the Ricci scalar, 
such as /(R) = R + a n R n + a m R m (where the constants a,j have the dimensions of 
length 2 ^' 2 ) has been investigated in [S], [5] and [TU] as far as cosmological implications 
are concerned. No analytical solution for the weak-field limit of the corresponding field 
equations exists, not even for (physically- relevant) special choices of the exponents n 
and m. 
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3. Non-analytical /(R) models 

In this section, we briefly review the main features of an example of non-analytical /(R) 
model, namely 

/(R) = R + 7 R /3 , (5) 

where the parameter 7 has the dimension of length 2 ^~ 2 ; and the exponent (3 is 

dimensionless. The typical length-scale L 7 of the model can be worked out of the 7 

1 

parameter as L 7 =| 7 I 2 / 3 - 2 . The interest in this model is based on the fact that field 
equations are analytically solvable, as in |2] and [TT] . 

The weak-field limit of (j5| was studied in [2J. The weak- field limit of field equations 
can be solved consistently as a perturbation of the flat metric only for 2 < (3 < 3. In 
fact, for > 3, the correction would be > 0(h 2 ), for (3 < 2 the correction would be 
grater than the ordinary 0(h) terms for R ~ 0, while for (3 = 2 the analytical case 
would be recovered. 

Because of the properties of the modified Einstein equations and of the functional 
dependence of ^ on R, it is straightforward to understand that the solution of the 
pertinent Einstein equations in the weak-field limit will consist of two parts, namely the 
standard GR term plus a correction term. In particular, in [2] , the spherically-symmetric 
Einstein equations were solved, and the generalized gravitational potentials were found 
to consist of a Schwarzschild term (responsible for the Newtonian behaviour) plus a 
modification term, whose features depend crucially on the ^R 13 term. The parameter 
space of the model was constrained by imposing compatibility with Solar-System data. 
As a result, a wide range for the values of the parameter 7, hence of the characteristic 
length scale L 7 , was demonstrated to exist, and a lower bound was determined. 

The cosmological implementation of this same model has been addressed in [TT] . 
In particular, the possibility to recover both a radiation-dominated era and a matter- 
dominated one has been demonstrated. From these calculations, also an upper bound 
for L 7 has been found. Combining the two results, it is possible to determine the 
allowed values for the characteristic length scale L 7 as a function of (3 (see figure Q, 
which was obtained in [HJ). Furthermore, the gravitational- wave evolution on a RW 
background has been illustrated to be, for all practical purposes, the same as in GR, 
and any contribution to the inflationary GW background was shown to stay below the 
detection threshold of present and future interferometers. 

In the following, we will analyze the weak-field limit of the model as far as the 
presence of gravitational waves is concerned. 

4. Minkowskian limit 

We will now solve the field equations in the weak field limit, and describe the main 
features of gravitational waves. To do so, we will recall the fundamental gauges of 
GR, which is possible to impose in order to get physical insight on the solution of field 
equations. After this, we will briefly recall the known results for analytical models, for 
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Figure 1. The gray zone represents the allowed characteristic- length scales of 
the model /(R) = R + 7R^. L^ Im comes from the Solar-System data and 
jMax comes f r0 m the study of the cosmological implications. See 



which a comparison is useful. We will then find the solutions of the field equations in the 
non-analytical case, and discuss them both on the basis of their physical interpretation 
and of their mathematical well-posed-ness. 

In the weak- field limit, we can evaluate the Ricci tensor, R^, and the Ricci scalar, R, by 
retaining terms that are first order in h^ u only. Thus, we get the expressions of the Ricci 
tensor R^ and of the Ricci scalar as a function of the perturbation h^, respectively: 

Ray [ E- l/lui/ ~\~ (h i/,aa h n,ua h ct.ujy)] , , . 

2 ' (6) 

R = -Dh + h a \ av , 

where □ = rj^d^du is the flat-space D'Alembertian operator [|] 
The Einstein equations ^ for our model ^ read: 

*V - Is^uR - 7/3V M V,R /3 - 1 + T/^g^g^VpV.R^ 1 = 0, 

3 7 /3g^V M V,R /3 - 1 - R + 7 (/3 - 2)R /3 , 
and, in the weak-field limit(§J they become: 

R»u - \%uR - iPV^vR?- 1 + wPUB?- 1 = 0, (8a) 

3-fPDRP- 1 = R, (8b) 

\ In the present paper, the definition of the flat-space D'Alembertian operator differs from that given 
in [3], where □ = —-q^d^d^. 

§ We remark that the weak field limit of the differential operator g^V^Vu (which appears in ^ acting 
on (the weak-field limit of) an object X such that O(h) < X < 0(h 2 ) is the flat-space D'Alembertian 



operator □ = rj^ v d^d v (as in (8a|-(86|) 
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This way, it is possible to rewrite Einstein equations (Sa)-(Sb) for the physical unknowns 
h^ v and h through he weak-field expansion (|6]). 

Analysis of the allowed gauges The features of the physical unknowns h^ v is at the 
basis of the determination of the possible gauges that can be imposed. 
We recall that the 'de Donder' gauge allows one to reduce the expression of the Ricci 
tensor R^ and that of the Ricci scalar R to 

1 2 (9) 

R = — nh 

2 

where □ is the flat-space D'Alembertian operator, by imposing g^T 01 ^ = 0, i.e. 
W 1 v „ = \h x X u . Because of the specific expression of field equations in General 
Relativity, this gauge is also known as the 'harmonic gauge' because Einstein equations 
in vacuum are solved by harmonic functions for h^. 

As well-known, in modified theories of gravity, the solution of field equations is not 
given only by harmonic functions. Nevertheless, because the de Donder gauge is based 
on the discussion of the degrees of freedom, we can apply the same reasoning to the 
modified field equations. 

Because of the linearity properties of the D'Alembertian operator □, it is always possible 
to split the perturbations into two parts, i.e. 

hp, = h<$ + h$, (10) 

where h ®} is a harmonic function, i.e. Oh$v = 0, while hpv is not, i.e. Oh^u 0. An 
analogous discussion can be made for h = rj^h^ = + h^K 

We remark that the non-trivial part of field equations (8a)-(8b) contain and h 
only through the expansion ([9]), i.e. they account only for the non-harmonic functions 
h$ and fcM. 

Furthermore, the sum of two harmonic functions is still a harmonic function, by which 
the harmonic part h$ can be suitably defined. Accordingly, because of the explicit 
expression of the field equations (8 a)- (8b) as functions of h^ v and h through ([9]), it is 
physically equivalent (apart from the order of the corresponding equations) to discuss 
the field equations for the physical unknowns and h or for the geometrical objects 
R^ v and R, defined as in 

The 'transverse-traceless' (TT) gauge can be imposed on h$ by using the four 
degrees of freedom left by the infinitesimal transformation of the coordinates, which is 
performed when fixing the de Donder gauge. Thus we get h^} TT . 



Analytical case We briefly discuss here the main results of the Minkowskian limit of 
analytical f(R) models. 
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The weak-field limit of field equations for the unknowns R^ and R read 
Rnu = qV^R + 2a2-R,^ , 
OR = jr^R, 



which are formally solved for the physical unknowns h 
taken into account) as 



(i) 



R 



Qo ~ Q 2 



rj^R + 2a 2 q ll q 1/ R 



fonce the solution for R is 



:i2a) 



[12b) 



A being an integration constant with the dimensions of length' 2 . 

The solutions (12a) and (126) of field equations imply the presence of a harmonic part 
h$ and a modified part h$ consisting of a massive mode with mass m, such 



that g 2 . — q 2 = m 2 , with several possible polarizations. It is interesting to remark, for 
future purposes, that the solution (12a) consists of two different parts: one containing 
the flat-space tensor i]^ u , and one containing the object q^q v . In particular, the latter 
is interpreted as a (constant) polarization tensor. Furthermore, this is obtained, in the 
solution of field equations, from the derivatives of the Ricci scalar (126). The functional 
dependence of the Ricci scalar (12b) on x M such that R = R(q^x^) is the only one that 
allows one to obtain a constant polarization tensor q^q v from the derivatives of R. 
Despite solution (12a) is obtained directly from field equations, in the literature different 
assumptions have sometimes been made, such as those in [12] and [15] . A mathematical 
classification of the possible (extra) polarizations that can appear in the weak-field limit 
of f(R) models, according to the features of the solution of the geometrical objects R^ u 
and R, has been performed in [13]. The behaviour of Weyl scalars in GR in the case 
of a binary black hole inspiral and merger has recently been numerically investigated in 
[15J. The comparison with the case of f(R) models could be a fruitful tool to further 
test the viability of modified theories of gravity in these scenarios. 

It is straightforward to remark that the scalar degree of freedom introduced by 
f(R) models allows one to rewrite the trace equation for the scalar degree of freedom, 
say ^, as 

□tf = m 2 ^ (13) 

where one identifies \I> = R, and the term m 2 \I/, with m 2 = 01/(602) the effective mass 
of the scalar field, as the linearization of a potential V(^/) (ruling the dynamics of 1 J r ) 
around its minimum, which corresponds to R ~ 0. The fact that an f(R) is analytical, 
i.e. its Taylor expansion holds in the neighborhood of R ~ 0, is therefore equivalent 
to linearizing the dynamics of the scalar degree of freedom $ on a fixed background, 
where the potential V(^) does admit a Taylor expansion in the neighborhood of such a 
minimum. 



This way, the choice of (126) accounts for a term describing a propagating massive field 
in (pah. 
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4-1. Non- analytical weak field limit 



We can now solve the weak-field limit of the Einstein equations for our model, (80 86). 
The trace equation (86) is solved by 

R=(6£)t&, (14) 

where the variable £ has the dimension of a length, the constant S has the dimension of 
length 1 '^ 5 , and read 



5 



(15) 



C*2 being an integration constant with the dimensions of a length, the scalar q^x^ = 
qo ct — q x x — q y y — q z z being defined with four integration constants q^. 
Substituting the solution (14) of the trace equation (8b) into the remaining equations 
(8a), and by taking into account (|9j), we obtain weak- field limit of field equations for 
the physical unknowns h^ u and h in the de Donder gauge, 

□/»$ = c^e a, 

nh {1) = 

which is straightforward integrated as 



(16) 



2(/3-l) 

/,(!) _ Q £p=2- 



where we have defined the objects Cn V and Q^u as 



a 



/IV 



l 2 n os .ze=}i2(p-l)p 
-^g^- 2 -2^7/35 k _ 2)2 , 



(17) 



;is) 



n 



(/3-2) 2 



6(/3-l)(g 2 - 9 2) 



1_\ P- 2 

2 7 



{ql 



P-2 f 1} P-2 



,2\ /3-2 



j^A f (/3-2) 2 
3(/3-l) 



(3-2 



(19) 



and the contractions C = rf^C^, Q = rj^ft^. 

The complete expression of ft,^ and h, defined in (10), is obtained by considering 



that h\ai and obey, by construction, the field equations 

□/$) = 0, 

Oh® = 0. 



(20) 
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Discussion of the solution of the modified Einstein equations We have solved our field 
equations by first solving the trace equations, and then substituting it in the remaining 
equations. In this paragraph, we discuss the features of the solutions of the field 
equations, from the points of view of their physical interpretation and their mathematical 
properties. 

It is possible to analyze the physical meaning of the trace equation as far as the 
presence of a scalar degree of freedom in modified-gravity theories is concerned. In fact, 
(86) rewrites as 



(21) 



with R = $3=1, as suggested by the presence of a scalar degree of freedom. Equation 
(21) is an inhomogeneous wave equation for the field $. As a wave equation, we are 
suggested to look for a solution of the kind $(t, x, y,z,) = $(£), with £ defined in the 
above. Within this point of view, we then interpret its solution, according to its Fourier 
transform, as a linear superposition of non-linear functions of the wave-packets given by 
the massive modes. 

We can infer that a function /(R) which is not analytical in the neighborhood of 
R ~ implies a scalar degree of freedom ruled by a potential which is not analytical in 
the neighborhood of the value R ~ characterizing the weak- field limit. 

So far, the solution of the trace equation (86) is not in principle unique. 
Nevertheless, we see that, for the choice R = R(£), (86) rewrites as an Emden- Fowler 
equation for the variable w = i?' 3-1 , i.e. [16] : 



d 2 i 



(22) 



Emden- Fowler equations admit the general power-law solution ([14]), and, for certain 
values of the parameters, they may also admit a particular (parametric) solution. For 
the functional dependence of ( 22 ) on the parameters, no particular (parametric) solution 
is known to exist, and (14) is therefore unique ([15]. [T7]). 

We can also go the other way round and consider that the trace equation (86) can be 
restated as a function of the unknown w = R 13-1 , for which [16] the "travelling wave 
solution in implicit form" reads 



Ci 



37/3 



wP^dw 



dw = q ct - q x x - q y y + q z z + C 2 



(23) 



For a generic choice of the integration constants C\ and C 2 , (23) is formally solved as 



i+ 



3 7 /3 2 C 1 (<^-q^) 



-rdw 



w 2 Fi 



2/3-1 . 
' 



2Q8-1) 



3 7/ 9 2 C 1 ( 9 2-g2 ) 



(24) 



For some special values 



where 2 F\ denotes the Gauss hypergeometric function 
of its arguments, the 2 F\ function can be expressed in terms of elementary functions, 
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but it is easily checked that, due to the range of the /3 values (2 < /3 < 3), in our case 



it is not [18]. Nevertheless, we notice that, by choosing C\ = 0, (23) admits the explicit 



solution (14). 



5. Interaction of modified gravitational waves with test particles 



We now use the geodesic deviation equation to study how the separation vector between 
two particles, A and B, changes because of the presence of a modified gravitational wave. 
In general, we have: 
D 2 



Dr 2 



;Sx" = R^ aM U p U a 5x\ 



(25) 



where 5x^ represents the separation vector between the two considered particles. 
A convenient coordinate system for analyzing the previous equation is the proper 
reference frame of one of the two particles, say the particle A. This frame has spatial 
origin x- 7 = attached to A's geodesic and time coordinate equal to A's proper time, 
hence x° = r on the geodesic given by x^ = 0; in addition, we assume that this frame is 
non-rotating (see [IB])- The coordinate system we have just defined is a local Lorentz 
frame all along A's geodesic, so that: 



ds 2 = dx° -5 ll dx i dx 1 + 0(\x j \ 2 )dx^dx a . 
Because of our choice of the reference frame, we also have 5x^ = x J b- 



(26) 



Moreover, at x 3 = 0, 
this geodesic, we find 



p vanish and so does dT fl s ■s/dr, thus, if we evaluate (25) along 



d 2 x j B 



R 3 



r k 



(27) 



dr 2 " k 

The curvature tensor is gauge invariant in the linearized theory, so we can write it by 



making use of the perturbation to the flat metric, given by h 
h^} TT and h$ defined in the above. As a consequence, we obtain: 



[IV 



h$ TT + h 



(1) 



with 



R 



k 



2 V ikfiO ik,00 



h {1 ±~ - h>tL^ -h.. 

OO.ifc iOfik Ok,iO 



,(1) 



,(1) 



(2* 



In order to simplify the notation, we now drop all the symbols 



which appear in 



(]27|) and (|28|). Eventually, we get: 

1 



d 2 x B 
~d^ 2 ~ 



, (0)TT , (1) 
u ik,00 ' u ikfiO 



+ h W 



h {1) -h W 
"-iOfik u 0k.i0 



(29) 



A key observation for solving (29) arise from the fact that we are treating a weak 



gravitational field, so we can also apply the slow motion approximation for the particles, 
so that: 



U ( 



U' 1 ~ 0, 



1,2,3. 



(30) 



Because of the relation (30), we can make the following assumption: 



X 



X T 



A' 



(31) 



Non analytical f(R) theory of gravity 



12 



where the A*'s are constant quantities. The scalar £ rewrites £ = (qod + £), where we 
have defined the constant S as S = (C2 — g^A 1 — q y A 2 — q z A 3 ). 

Furthermore, we assume that, initially, the particles are at rest relative to each other, so 
and we consider the particle relative position displacement, induced by 



that x^ 



x 



B(0) 



the perturbation h^, as a small perturbation with respect to the initial position x 
Due to this fact, we can set: 



B(oy 



x 



B 



X 



B(0) 



+ X 



B(l)l 



(32) 



where x J B ^ is a constant vector, and the index j refers to spatial coordinates only. 



Collecting all the ingredients together, solving ((29]), (j32J) rewrites 

J 



x 



B 



X 



B(0) 



ik 



X 



B(Q) 



+ 



00 ■ 



/A: 



.a 



Ok 



where: fc = 1, 2, 3, thus & = -<f = -q z ) and g = g' 



2tf-l) 



(33) 



B(0)" 



Geodesic displacements To study the features of the modified gravitational waves in 
our model, we choose particular initial conditions for the wave, such that the effects 
of the modified gravitational wave are easily pointed out. Furthermore, we will choose 
the most appropriate numerical values of the parameters, such that the constraints on 
the modified part of the gravitational wave are the strictest. We will also illustrate the 
effects of the extra polarization modes of the gravitational wave. 

To do so, in particular, we assume the h^v part as propagating along the x direction, 
in the A + polarization, and the h$ part as propagating along the y direction. The 
relative displacement of the particle B under the influence of the gravitational wave, 
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(34c) 



where the terms containing £ come from h^ 1; 



/IV ■ 



From the previous relations, it is possible to check that we are allowed to define the 
velocity v of the GW modified part h)S, as v 2 = c 2 q\jq 2 y . In order to obtain a numerical 
estimation for v, we now assume to reveal GWs today, i.e. at t = t u = 14 x 10 9 years, 
and at y = x 2 = 0. Moreover, we note that £ = q^x^ + C2, where C2 is an arbitrary 
constant, then, we also make the assumption that, for our choices, we can neglect C2, 
or set C2 = 0. 
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Since there are not direct evidences of the gravitational waves, like, for instance, 
observations by means of ground-based interferometers, we can also assume that the 
modified part of the GWs in our f(R) gravity cannot produce relative displacements 
larger than those produced by the ordinary gravitational waves. If we now consider 
A + ~ 5 • 10 -23 for the ordinary part h$ TT jjj we can associate with this term relative 
displacements of the order of 10 -23 . Thus, setting /3 = 2 + n, we obtain the following 



expressions v 2 from (346) and from (34a) and (34c) 

r 2 

2 ° 

" 2 = i+ 5 



(35 a) 



XQ — 23n 



Vl,3 



1 + 



(356) 



1 



For the picture considered above, after fixing an allowed value for L 7 , (35a) and (356) 



depend only on f3, thus it is possible to determine numerical estimations for both v-i and 
Vxfi- Hereafter, we will take into account an intermediate value for the characteristic 
length scale of our model, i.e. L 1 = 38pc. It is worth noting that, in the limit — > 2, we 
obtain that both v 2 and t> lj3 tend to c, while, for /3 — > 3, we get v 2 ~ 6.197 x 10~ 21 km/s 
and ui f 3 ~ 4.3 x 10~ 8 A;m/s. For the discussion in the abaove, (35a) and (356) give us the 



maximum value the velocity v can have, so that we do not obtain relative displacements, 
due to hp}, larger than those expected from the ordinary gravitational waves. Because of 
the different dependence on n and because of the presence of the factor c 2 ™, it is possible 



to check that (35a) requires a smaller value of the velocity with respect to (356). As a 



consequence, we have to assume as a good estimation for v, the values obtained from 



(35a), indeed it imposes more restrictive constraints. 



Modified polarization We now want to show that the modified gravitational waves, 
described by the two terms h^v and hpi, change the polarizations of the standard 
gravitational waves. Here, we will keep on assuming that the term hp°} TT describes an 
ordinary gravitational wave propagating along the x-axis with plus polarization only and 
the modified part of the gravitational waves, indicated by h^}, is charactered through 



the quadrivector g M = (q ,0,q , 0). With these assumptions, (34a), (346) and (34c) are 



still valid; furthermore, if we consider the arbitrary constant C 2 as negligible or set it 
to zero, we have: 



f = cq (t-y/v). 
According to (9) of [19], where we are using G — 6.67 x lO 11 ™ 3 ^^ - 



(36) 

and c = 3 x 10 8 m/s, 



it is possible to obtain an amplitude of the standard GWs of about 10 , if we use the following 
approximated values for a binary star system: Object masses = M Q , Orbital radius = 20 x 10 3 m, 
Orbital frequency = 400Hz, Source distance = lOOMpc. 
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By making use of the relations (34a), (346), (34c) and (36), it is possible to verify that 
we can define the argument of the GW modified part as the dimensionless quantity 
(c/L 7 )(t — y/v). We now observe that, in our calculation, we are going to use the 
previous velocity estimations. As a consequence, we point out that we estimated the 
velocities of the GW modified part, by imposing that, at y = and t u = 14 x 10 9 years, 
the relative displacements due to the term h^} were of the same order of the relative 

(0)TT 

displacements associated with h^v ■ Since we want to show the effects introduced by 
the term h^} on the ordinary gravitational wave polarizations, we need to consider only 
those situations where the observable effects due to the term h$ are at least comparable 
with the observable effects due to the term h^} TT . Hence, we will take as reference value 
for the argument (c/L 7 )(t — y/v) of the GW modified part the one we have used when 
deriving the velocity estimations, that is: 

c f t _V\ =^1 = 1.16184x108, (37) 



L 7 V v ) L 

where we keep on using L 1 = 38pc. As far as the relative displacements due to the term 
^(o)tt concerned, for the particular conditions we are considering, we will use the 
following definition: 

,(0)TT ,(0)TT 

H—H"^ (<"('-;))■ (38) 

with h + = 10 -23 . To visualize the effects introduced by the term hf^o on the GW 
polarization, we have taken into account a system of particles initially arranged on a 
circumference, on the plane x = 0. Here, we insert the plots for this case (see figure ^ 
and figure (l3j)), where we have amplified the effects due to the modified GWs, for the 
sake of an effective comparison. 

In figures ^ and (|3]), as far as the argument of the modified GW is concerned, we 
have adopted the same values for both h^} TT and h$. 

Our results tell us that, in the case of /(R) = R + 7R^ theory, we obtain polarization 
ellipses which change with time, similarly to the General Relativity case. However, 
differently from General Relativity, in our /(R) gravity theory, the polarization ellipses 
keep changing and do not assume periodically the same shape. We also note that the 
larger the argument values, the greater the effects of the GW modified part. This last 
result comes from the fact that the amplitude of the modified GW term associated with 
h$ increases when the argument increases. 

Moreover, from (34a)-(34c), we can see that, in general, modified GWs have effects 
along all the three spatial directions, then GWs, in /(R) = R + 7R^ gravity, are not 
transverse. 



6. Concluding remarks 



In this paper, we have analyzed the features of the model /(R) = R + 7R /3 , whose 
Taylor expansion in the weak-field limit R ~ does not hold. 
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Figure 2. Left column: polarization ellipses in General Relativity. Right 
column: polarization ellipses for f(R) = R + jR^ gravity theory when 
/3 = 2.0000001 and L 7 = 38pc. For the sake of clearness, we note that we 
have used the same argument 3l7r + (H — 100) < p =< 33ir + (H — 100) for both 
the two terms of the modified GWs and GWs in General Relativity. 



In Section |2j we have pointed out the main features of modified f(R) gravity, focusing 
our attention on the characteristics of the weak-field limit, as far as analytical models 
as well as non-analytical ones are concerned. 

In Section [3j we have reviewed in some detail the properties of the non-analytical model 
/(R) = R + 7R/ 3 , which have already been investigated in previous works, i.e. [2] for the 
weak-field limit in the spherically-symmetric static case, and |TT] for the cosmological 
implementation. 

In Section |4j we have found the analytical solution of the weak-field limit of field 
equations as far as the presence of gravitational waves is concerned. To do so, we 
have briefly recalled the properties of some useful gauges in GR, and their effects on our 
modified field equations. After this, we have summarized the main results found in the 
case of an /(R) model whose Taylor expansion in the weak- field limit around the value 
R ~ holds. The explicit solution of the field equations accounting for gravitational 
waves has then been given, and its features have been investigated and commented as 
far as their physical interpretation is concerned, as well as its geometrical structure and 
its mathematical well-posed-ness. 

In section [5j we have studied the interaction of the gravitational waves described by the 
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General Relativity 
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Figure 3. Left column: polarization ellipses in General Relativity. Right 
column: polarization ellipses for f(R) = R + jR 13 gravity theory when 
(3 = 2.0000001 and L 7 = 38pc. For the sake of clearness, we note that we 
have used the same argument 33-7T + (S — 100) < p = pgr < 357r + (S — 100) 
for both the two terms of the modified GWs and GWs in General Relativity. 

/(R) = R + 7R^ model with test particles. More precisely, we have solved the geodesic 
deviation equation and have imposed some constraints on the velocity of the modified 
gravitational wave. Furthermore, we have discussed the polarizations which are present 
in this model. 

Thus, we have shown how the particular modification of General Relativity we addressed 
here is able to provide a very peculiar trace of its presence by means of the emergence of 
a non-linear and non-negligible correction to the linear theory of the gravitational waves 
propagation. Such a specific morphology of the modified spacetime ripples, more than in 
their non-transverse and non-traceless behaviour (present in other modified gravitational 
wave paradigms) , consists in the details of their polarization and propagation features. 
In fact, our study of the polarization tensor, together with the increasing power-law 
behaviour of the modified term, induce a well-determined picture of the polarization 
ellipses, which seems the very smoking gun of the modified power-law Lagrangian we 
added to the standar Einstein-Hilbert shape. This striking morphology, well-grounded 
on the amplification that any signal of this sort acquires in a given space point, make the 
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issue presented above of some impact also for the expectation of the future detection 
of gravitational signals, even because the proposed non-linear spacetime ripples are 
specifically traced in the stochastic sea of all the other astrophysical signal, at the 
lowest order, properly accounted by General Relativity. 
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